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hazard further conjecture as to the reason o£ this difference of 
behaviour of the instrument in the two elements we might 
suppose more definitely that the level error is caused by a 
warming of the eastern pier more rapidly than the western, 
whose conductivity is not so great; so that it expands vertically 
and affects the level, while at the same time the lateral expansion, 
being symmetrical about the central vertical plane which passes 
nearly through the pivots, would not affect the azimuth. The 
piers, or that part of them at least which supports the pivots, 
not being very thick, the lagging due to gradual conduction will 
not be very large, and we have the effect of difference of conduc¬ 
tivity nearly free from any counteraction. 

On the other hand, we must refer the azimuthal variation to 
changes of temperature probably at some depth below the 
surface of the soil, which do not happen to affect the level 
appreciably. 

I have examined the variations of level and azimuth at other 
observatories for a few years to see if there is any similarity 
between these fluctuations for different instruments. It is 
somewhat difficult, without expending more time than I can 
well spare at present, to disentangle their annual variations 
from changes, secular and irregular, and generally larger than 
those noticed at Greenwich. I look forward with interest to a 
discussion by Mr. Finlay of the errors of the Cape Transit 
Circle promised in the last volume of “Cape Observations” 
(1879-81). Roughly speaking, there does not seem to be much 
similarity between the fluctuations of position of different instru¬ 
ments ; for instance, the difference of epoch, which has been 
chiefly considered above, is very variable. And indeed Mr. Ellis 
pointed out in his paper above referred to how largely the 
variations of the present Transit Circle differ from those of 
Troughton’s Transit, which occupied nearly the same site. As 
I have already said, however, in suggesting the above explana¬ 
tions of the changes in position of the Greenwich Transit Circle, 
we have at our disposal two causes which, combined in varying 
relative proportion, are capable of explaining fluctuations follow¬ 
ing or anticipating the temperature by very different periods of 
time. 


On the Formulae for Computing the Apparent Positions of a Satellite , 
and for Correcting the Assumed Elements of its Orbit. By A. 

Marth. 

The methods of computation connected with the investiga¬ 
tion of a satellite’s orbit are comparatively simple and convenient, 
if they are duly selected to suit directly the coordinates furnished 
by the observations. If polar-coordinates have been observed, 
there exists apparently no good reason for not employing polar- 
coordinates also in the corresponding computations. The present 
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paper contains the formulse, which seem to me the most advan¬ 
tageous and handy for computing the apparent positions and 
the corrections of the assumed elements in the chief cases which 
come under consideration. 

Referring the plane of the satellite’s orbit to the plane parallel 
to the terrestrial equator, let N denote the longitude or Right 
Ascension of the ascending node ; J the inclination ; Q the angle 
between the line of apsides and line of nodes or the departure of 
the lower apside from the node ; e = sin <f> the eccentricity : p,, e, v 
the mean, eccentric, and true anomaly of the satellite; u = Q4-v 
its true longitude in the orbit reckoned from the node ; u Q = Q+ 
its mean longitude ; r the radius vector ; a the major semi-axis of 
the orbit or the satellite’s linear mean distance from the planet’s 

centre; ^ = p, so that p is the radius vector expressed in parts of 

CL 

the major semi-axis of the orbit. Further, let A, D denote the 
apparent Right Ascension and Declination of the planet’s centre ; 
a, 8 those of the satellite ; A the linear distance of the observer 
from the planet. If the origin of the coordinates is placed in the 
satellite instead of in the centre of the planet, the strict formulae 
become somewhat simplified, but for the purpose a and 8 must be 
known. In case approximately correct values of a — A and 
8 —D are available, like those of the three outer satellites of 
Saturn , published for some years past in the Monthly Notices , or 
if they can be derived directly from the observations, it will be 
worth while to take advantage of the favourable circumstance. 

As most of the observed positions of satellites consist of 
measures of position-angles p and distances s, or are given in 
polar-coordinates, the formulae for their proper computation may 
first be taken into consideration. 

If the Earth’s position, as seen from the satellite, is referred 
to the plane of the satellite’s orbit by the longitude IT + i8o° and 
latitude B, and if P is the position-angle of the pole of the orbit 
(in Right Ascension N—90° and Declination 90 0 —J) at the geo¬ 
centric place of the satellite, the values of U, B, P are found from 
those of N, J, a, S by means of the equations : 

cos B sin U = cos 8 sin (a —N) cos J + sin 8 sin J 
cos B cos U = cos 8 cos (a—N) 

sin B = cos 5 sin (a —N) sin J — sin $ cos J 
cosBsinP = — cos (a—N) sin J 

cos B cos P = sin $ sin (a—N) sin J + cos 8 cos J 
or by equivalent formulae. 

The satellite’s orbital longitude u=Q-\-v, and also p being 
deduced from the mean anomaly p, either by the usual formulae: 

= p v> = 57*296 

• • J1 J am 1° 

p sm v = sin € cos </> 1 

p cos r = eos e — * 
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or by the series for the equation of the centre v—pc and log p , 
the values of the position-angle p and of the distance s corre¬ 
sponding to u are found by means of the equations 

sin <r sin (p —P) = sin (u — U) 
sin crcos(p—P) = cos (zc — U)sinB 

. a 

s m s = ~ . psm cr 
A r 

in which o- is the arc on the sphere between the geocentric and 
planetocentric place of the satellite, or, without the introduction 
of O', by 

sin 5 sin p — P) = !L p sin [ii — U) 
sin s cos (p — P) =£ p'cos(w—U)sinB 

The position-angle p thus computed refers to the Declination-circle 
passing through the satellite. In ease the observed position- 
angle refers to the Declination-circle passing through the point 
midway between the planet’s centre and the satellite, a correction 
= -— ^-ssinp tan 8 must, when a strict comparison is required, 
be added to the computed p. 

If A, D are substituted for a, 8 in the formulas for U, B, P, so 
that their ensuing values refer to the centre of the planet 
instead of the satellite, the correction of the computed p will be 
= + \ s sin p tan D, the computed <r (which must be taken 
between 90° and 180°, when cos (u— IT) is negative) represents 
the arc between the geocentric place of the planet and the planeto¬ 
centric place of the satellite, and the value of s must be found 
indirectly from the equation 

sin s=—p sin(cr — s), 

A 

if the employment of the equations 

Aj sin s sin(^ —P ) = ctp sin(w — U) 

A 1 sin s cos [p — P) = ap cos (u — U) sin B 
A, cos s — ap cos (u — U) cosB-r A. 

in which Aj denotes the distance of the satellite from the earth, is 
to be avoided. 

The equations of condition for correcting the assumed ele¬ 
ments of the orbit gain in clearness and simplicity by the intro¬ 
duction of 

sin o’ sin r = sin B 
sin cr cos r = cos B sin (u — V) 
cos 0 - = cos B cos(^ — U) 

which serve, at the same time, as a partial control upon the 
equations for finding p : 
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sin o sin (p — P) = sin (u — U) 
sin o cos ( p — P) = cos (u — U) sin B 
a 

sm 5 = — . p sm o 

A 

The differentiation of these equations gives : 


sin o'dp = sin r . d {u — U) — cos r cos (u — U) dB 4- sin odF 
do = cost . d(u-~ U) + sinrcos(w — U)dQ 3 

cos BtZU = — cos U sin B dJ % — cos P cos &ZN 
cosIMP = — cosTMJ 4 - sin U sin JtZN 

d B= +sinU^J — cos U sin JcZN. 


Hence, after some reductions : 

sin odp — sin rdu — cos r sin udJ 4- (cos o cos p eos 8 — sin o sin 8)tM 
do — cos rdu 4 - sin r sin ud J 4 - sin p cos SdN 

p sin o . «° — = p cos odo 4- sin ou° dp + p sin ooo° . —. 

s u 

The differential expressions of u, the satellite’s true longitude, 
with respect to its mean longitude u Q and Q and <£, are : 

pdu = ( ML c tdu 0 4- (p - 4- (1 4 - p sec ? 0 ) sin v . cos <pd(p 

P V p ) 

— (1 4- e cos v ) sec <pdu 0 — {2 cos v } ten <pdQ 4- {2 sin t} cos <pd<p 


in which the coefficients in {} brackets denote the following 
equivalents, among which the most convenient, according to 
facilities at hand, may be selected : 


[2 cos v } = cos v 4- cos e cos <p 4-tan 

_ (2 4- e cos v) cos v 4- e 4 - tan cos 2 <p 
I + e CoS v 

e cos (f> + tan 


~ 2 - 


e eos v 
1 4 - e cos 


iVj 


cos v 4- 


1 4- eccs v 


= (1 4 -p sec. 2 4 >) cos v 4 - p(tan <p 4- tan %cp sec 2 <|>) 
{2 sin t } = sin v 4- sin € sec <p 
2 4- e cos v 


_ sm v 


1 + e cos v 

( eeosv 


= I 2 


i 4- e cos to 
= (l +psec 2 cp ) sin v. 


sm v 


The substitution of these expressions and of 

o >°dp — tan cf> sin vdu Q — sin v . tan cpdQ, — cos v . cos <pd<$> 

in the equations for sin crdp and ds leads to the equations be- 
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tween the variations of the computed coordinates and the varia¬ 
tions of the assumed elements 


p sin <rBp = sin r . . Su Q 

P 

— sin r. {2 cos r] w°(tan (p + 5 -) sin 5 Q 

+ sinr . {2 sin t’}ar[(taii <p - 5 ?) cos 5 Q —tan <f>] 

— p cos t sin 

4 - p(cos <t cosp ccs 5 - sin cr sin 5 ) 5 X 


P sin cr. ct>°~ = (cos a cos r . 1^— 4 - tan <p sin cr sin r)5z/ 
5 p 


— (cos <r cost{ 2 cos v] - sin <r sin r)co°(tan $ + 8 e)sin SQ 
+ (cos <r cost{ 2 sin v) —sin a cos c)& 0 [(tan$ + &2)cosSQ — tan<£] 
+ p cos <r sin r sin zz. 5 J 
+ p cos o' sin p cos 5 . 5 X 

+ p sm (T . w°— . 

a 


In order to guard against some possible oversight or error 
in the interpretation of the variations connected with the eccen¬ 
tricity if they are written in the form tan 6 0 Q and cos <£ 8 cf>, I 
have substituted in these equations for 


and for 


tan <pBQ .w°(tan <p + Be) sin SQ 

cos (p$(p .... «°[(tan <p + Be) cos 5 Q — tan < p] 


expressions which represent the true meaning of these variations, 
and will lead to right deductions of the corrections of Q and e in 
all cases. If not thus interpreted the expressions tan 8 Q and 
cos </>$</> or Se may obviously be misleading, and may give very 
incorrect results whenever the assumed values of Q and e 
require more than small corrections. 

If, instead of 8 Q and 8 e, the variations of suitable functions of 
Q and e are to be introduced, the equations for p du and dp suggest 
the advantage of selecting for the purpose not e sin Q and e cos Q, 
but E sin Q and E cos Q, where E is such a function of e or <£ that 


tan <pd E = E de = E cos (pd<p. 

Putting accordingly 


E = e . sec 2 \<p . c ~ 2 sin2 ^ = 2 tan ±<p. c - 2 sin= (c basis of nat. log.) 


the differential coefficients of u and p with respect to E sin Q and 
E cos Q will be 


pdu — 


(2 + e cos v) cos u+(e cos $ + tan l 4 >) cos Q 
1 4- e cos v 


EsinQ) 

E 


+ 


(2 + eeosv) sin u + (e cos <£> + tani<£) sin Q 


1 + e cos v 


tan (p 

E 


d(E cos Q), 


dp-. . . . — sin zz . . d( E sin Q) — ccs u , —. d\ E cos Q). 
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Hence the eccentricity terms in the equations of condition for Sp 
and S s become 

* . (2 + e cos v) cos u + (e cos d> 4- tan ) cos Q 0 

p sm crop = .,.-8111 r . i- 1 -_ . o>° 


i + e cos y 


151?. 5(E sin Q) 


+ sin -r ( 2 + e cos v ) si n u + (e cos <p 4- tan ±(f>) sin Q /<Q 

1+ € COS V 


85 

p sm <r — 


l^.S(EcosQ) 

, # — (cos a cos r + e CQS v '1 cos u + ( e cos ft + taD Ml cos Q 


1 + e cos 

+ sin <r sin w) . . 8(E sin Q) 

E 


+ (COS <7 COS T . ( 2 + * C0S ^ sin M + ( g C0S !t + tan M sin Q 

I 4- e cos -y 


— sin a cos w). . 5(E cos Q). 

E 

Though in this way the corrections of the assumed values of 
Q and e might be properly found and mistaken deductions 
avoided, the way would be needlessly circuitous ; for, by substi¬ 
tuting in the foregoing expressions for 8(E sin Q) and 8(E cos Q) 
their values 

(E 4 5E) sin (Q, 4 5Q) — E sin Q, and 
(E 4 SE) cos (Q 4 8Q) — E cos Q, 


and by eliminating Q, the expressions are obtained which are 
already given in the equations of condition, and which, while of 
equal strictness, are obviously preferable on account of their 
greater simplicity and directness. Moreover, they allow the 
motion of the apses, if approximately known, to be easily taken 
into account. 

The equations of condition have yet to be so adjusted that 
the differences between the observed and the computed co-ordi¬ 
nates are measured by the same unit. It will depend on the 
character of the measured position-angles and distances whether 
it is preferable to express the differences in seconds of arc or in 
units corresponding to an arc of one degree on a circle of suitable 
radius, and also how far the ensuing rules are better modified or 
altered. 


Ifsin<z 0 =-^-, where A 0 is some conveniently fixed value of A 
(not differing greatly from the average of the distances A at 
which observations are taken), and if r, the angular distances 
s of the satellite will be 


s = a 0 v . p sin cr. 
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Making 

573 

the adopted unit, the equations of conditions will he 

-5j)=+j'smr. — . ou 0 

«0 P 

— v sin t . {2 cos v} . w° (tan 0 + Be) sin SQ 

+ v sin r . {2 sin y} . w° [(tan <p + 5e) cos 5Q — tan 0 

— vp cos r . sin n . 5J 

+ vp (cos <r cosp cos 5 — sin <r sin 5) . 5N 

5s , cos 0 , . . . v * 

_ = 4- j/(cos <r cos r . - - 4 - tan <p sin cr sin v) , ou 0 

K P 

— i/(cos cr cos r . {2 cos y] + sin <r sin v) . o>° (tan 0 + 5e) sin 8 Q, 

4 - v(cos <r cos t . {2 sin y}- — sin <r cos v) . a° [(tan 0 + Be) cos 8Q — tan 0] 

4 - yp cos o* sin t . sin u . 8J 
4 - vp cos <r sin p cos 8 . 8N 

af k‘ 

If the differences between the observed and computed coor¬ 
dinates are expressed in seconds of arc, the equations become 

sinT . C0S -~ . kBu q 

57 '3 ^ P 

~ v sin t . {2 cos v\ . a 0 (tan 0 4 - Be) sin 8Q, 

+ r sin t . {2 sin v] . a Q [(tan 0 4 - Be) cos 8Q — tan 0] 

— ypcosr . sin u . kSJ 

+ vp (cos (T cos p cos 8 — sin <r sin 8) . k8N 

5s ■■= 4 - v (cos o' cos r . 4 - tan cp sin cr sin v) . fr8w 0 

P 

— v (cos o' cos t . {2 cos y} 4 - sin o' sin v) . a 0 (tan cp ~ 5c) sin 5Q 

4 - v (cos o- cos r . {2 sin v } — sin cr cos v) . c? 0 [(tan 9 - Be) cos 5Q — tan 0] 
4 - vp cos <r sin r . sin u . tcBJ 
+ vp cos o* sin jo cos 8 . «5N 

4* — • 5c£ 0 . 

^0 

The formulae are here written out in full. Abbreviated 
notations for the quantities which are to be determined and for 
their coefficients will, of course, be substituted in actual com¬ 
putations. 

The preceding equations are valid for orbits of any assumed 
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ellipticity. For assumed circular elements they get considerably 
simplified, and become 

— 5p = + v sin r . 8 u 
a 0 

— v sin t . cos u . 2 <J> sin Q 
+ v sin t . sin u . 2<p cos Q 

— v cos t , sin zc, 8J" 

+ if (cos o' cosp cos 8 — sin o' sin 5) . 8N 


— = + V COS O' COS T . du 
K 


— v (cos o* cos r . cos u + sin <r sin u) . 2<p sin Q 
+ v (cos <t cos t . sin u— J sin o' cos u ) . 2<p cos Q 
+ v cos o' sin t . sin u . 8J 
+ r cos o- sinp cos 8 . 8N 


They must be multiplied by k if the differences are to be ex¬ 
pressed in seconds. 

The orbital longitudes u, U, Q employed in the preceding 
formulae are reckoned from the node N, and this reckoning is 
most suitable whenever the inclination J is considerable. But if 
J is only of moderate amount, as in the case of the orbits of 
Saturn's satellites, the determination of the nodal point is con¬ 
siderably uncertain, and this uncertainty rapidly increases with 
the decrease of J. By adding the longitude of N to the longi¬ 
tudes reckoned from N the uncertainty is evaded, and the 
longitudes in the orbit start from a properly fixed point. 

Putting u + 'N — l andU + ]Sr = L, the formulas for finding p 
and the sin and cos of o-and r require merely the substitution 
of Z—L for u — U : 

sin <r sin (p —P) = sin (Z —L) 
sin (T cos (p — P) = cos (l —L) sin B 

a . 

sin s = —. p sin <r or s = a 0 v . p sm <r 

A 

The differential expressions for sin crclp and dcr become by the 
substitution of dl —dN for du 

sin a dp = sin r . dl — cos t sin u . d J + (cos <r cosp cos 8 — sin o' sin 8 — sin 
dcr — cos t . dl + sin t sin u . dJ + (sin p cos 8 — cos 

But 

cos rr cos p cos 8 —sin a sin 8 = sin r cos J + cos r cos u sin J 
sinp cos 8 = cos r cos J — sin r cos u cos J. 

Hence 


sin <rsinr = sin B 

sin o' cos r— cos B sin (l— L) 

cos o' = cos B cos (l —L). 


sin adp = sin r . dl + [cos r cos ( l —N) — tan J.T sin t] sin J. dN —cos r sin (l — N)dJ 
d<r = cos t . dl — [sin r cos ( l —N) + tan \J cos r] sin J. cZN + sin r sin (l—~N)dJ 
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The equations of condition between the variations of the coordi¬ 
nates and the variations of functions of the elements become 
accordingly 



57-3 * 


— v sin r. {2 cos v} . a 0 ( tan (p + 5e) sin 5(Q, + N) 

4- v sin t . {2 sin v }. a 0 [(tan </> -t- 8<?) cos 8 (Q + N) — tan <f>] 



sinSN 


cos (t cos t . f tan (p sin o'sin . kS1 0 

— v (cos <r cost . {2 cost’} + sin <rsin v ). <z 0 (tan<£ + 8e)sin8(Q + N) 

4- v (cos <r cos t. {2 sin r} - sintr cos v ). a 0 [( tan (p + Se)cos 8(Q + N) — tan<^] 



+ —. 8& 0 
^A 


In order to guard against incorrect deductions of the cor¬ 
rections of the node and inclination, which might occur, if the 
terms sin J . SN and 3 J were kept without proper interpretation, 
I have substituted for 



and for 



expressions which, as in the analogous case of the eccentricity 
terms, represent the true meaning of these terms, and cannot 
mislead. By this simple substitution all the advantages for the 
correctness of the deductions are secured, which otherwise might 
be gained by the introduction of the variations of tan l J sin 1 ST 
and tan cos "N, but without their attending drawbacks, as 
mathematical readers may easily satisfy themselves. 

The quantities 

k 8/ 0 , <z 0 (tan<£ + 5e)sinS(Q + N), <x 0 [(tan<^ + 8e)cosS(Q -X)—tan$], . 



which are to be treated as the unknown quantities, will in the 
formulae for actual computation, of course, be represented by 
single letters. 

For assumed circular elements the equations of condition 
become 


c c 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at Mount Allison University on June 30, 2015 






1887MNRAS..47.. 333M 


34 2 


Mr. Mar Hi, Formulae for Computing 


xlyii. 6, 


57*3 


%> = + i/sinT. k5Z 


— j/sinT. cos l . 2K<p sin(Q, 4 - N) 

4- v sin r. sin l . 2k</> cos (Q, 4- N) 

( JJv 

sin J 4-—jsin 5N 

—x/cosrsin(Z —N).a 0 [^sin J + ^^cos5N — sin J] 

5s = 4- v cos cr cos t . /c5Z 0 

— v (cos <r cos t , cos l 4- J sin <r sin l) . 2 t«p sin (Q, + N) 

+ v (cos or cos r ♦ sin ^ sin or cos l ) . 2 K<f> cos (Q + N) 

( 5J\ 
sin J + j sin 5N 

( 5 J\ 

sin J 4- — \ cos 5N — sin J] 


+ _ . Sa 0 . 


In turning from the consideration of the proper treatment of 
measured polar coordinates to that of rectangular coordinates, 
it is obviously expedient to take first the special case of Saturn's 
satellites, when their rectangular coordinates parallel to the axes 
of the ring are either measured or estimated. In measuring the 
distances from the axes it is essential that the position-angles in 
which the measures are taken should be distinctly stated, so 
that if they differ from the computed position-angles of the axes 
which represent the fundamental plane, to which the orbits of the 
satellites are to be referred, the necessary corrections may be 
applied. 

If L, B, P have the same significance as before in reference to 
the assumed plane of the ring, the coordinates x", y n parallel to 
the axes of the ring are 

s sin (p — P) = x n = UqU . p sin (l~ L) 
s cos (p — V)= y" = a 0 u . p cos (l —10 sinB 
5#"= + v [cos(Z — L) + ecosM] . k sec</>5Z 0 

— r[cos (l — L)(cos ecos <p + tan + cos M] . a 0 (tan</> 4- 5e) sin 5 (Q, 4-N) 

4- v [cos (Z—L) sin e sec <}> + sinM] . a 0 [(tan <£4- 5e) cos 5 (Q, 4- IN') — tan <£] 

( 5 J \ 

sin J 4- __ J sin 5N 


4- —• . 5(X 0 . 
a n 


dy' — —v sin B [sin (l — L) — e sin M] . k sec <p8Z 0 

4 - v sinB [sin (I—L)(cos e cos (j> + tan^) — sinM].a 0 (tan <f> 4- 5<s) sin 5(Q 4- N) 
— v sin B [sin (Z — L) sin e sec (p + cos M]. a Q [(tan <p 4- 5^;) cos 5 (Q 4- N) — tan <p ] 

4 -^p[cosBcos(Z—N)-tan JJsinBsin(Z~L)] . a 0 ^sin J 4-^~^sin5N 

( 8T\ 

sin J 4- ) cos 5N — sin J] 

v" 

+ — . §« 0 . 
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For convenience M is put instead ofL — (Q-flN’). If the eccentric 
anomaly c is not available 

eosv + ecos<b +ta,n^<j> , 

— -————-must be put tor cos e cos <p -f tan 

I + e cos v 

and 

sin v j. 

- - for sin € sec <p. 

i+e cos v 

For assumed circular elements the equations of condition 
become 

8x" — + vcos (l — L) . k$1 

— v [cos (/ — L) cost + cos L] . a 0 c sin (Q + N) 

+ v [cos (l— L) sin l + sin L] . a Q e cos (Q fN) 

— v tan JJcos {l— L) . c 0 ^sin J + sin 5X 

x" * 

+ — . Sa 0 . 

a Q 

8y n — — ysinB . sin(Z —L) . k81 

+ j'sinB [sin {l — L) cos l — sinL] . tf 0 esin(Q, -f N) 

— v sin B [sin (l — L) sin l + cos L] . a Q e cos (Q, + 2s) 

— r [cosB cos (£—N)--tan-|J sin B sin (/—L)] . « 0 ^sin J sin SN 

+ vcosB sin (l— N) . a 0 [^sin J + cos 5N —sin J] 

a 0 

In case the x ordinates are not measured micrometrically, but 
estimated by comparison with certain definite points of the major 
axis of the ring*, the proper unit for x and for the semi-axis of the 
satellite’s orbit will obviously be the semi-diameter of Saturn’s 
equator. If (a) is the length of the semi-axis expressed in such 
units, the computed 


x~(a)p sin (l — L) 

and the equations of condition require merely the substitution of 
(u) for a Q v. 

Though the range of such estimations of x is limited, there is 
the great advantage that they can be made without micrometer, 
and even with a rather unsteady instrument, and that when 
taken with due care and circumspection at the right times, they 
yield in accuracy only to the best micrometrical measurements. 
Observations which fix the times of the conjunctions of the 
satellites with the limbs of the ball and the ends of the ring in 
the four quadrants are of great value, not only in investigations 

CC2 
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of the orbital longitudes of the satellites, but also in examining 
the dimensions of the ball and ring. 

If the observed positions of the satellite in reference to the 
planet’s centre consist of differences of Bight Ascension and 
Declination, or are given in equivalent rectangular coordinates, 
the frequently used method for their computation may be 
followed, which Bessel employs in his investigation of the orbit 
of Titan (“ Astr,-hTachr.,” Hos. 193--195). His auxiliary angles 
/, F, g , G, h, H, referred to the geocentric place of the planet 
A, D, may be found by means of the formulae * 


sin/sin E = — sin (A — N) 
sin/cos E = + cos (A — N) cos J 
cos/= — cos (A — N) sin J 
sin/sin E = — sin (A — N) sin J 
sin/cos E = 4 - cos J 


sin g sin (G — F)= —sin (D + E) 
sin *7 cos (G— F) = +cos(D + E) cos / 
cos g — + cos (D + E) sin / 
sin h sin (H —E)= +cos (D + E) 
sin Ji cos (H - E) = — sin (D + E) cos/ 
cos h = + sin (D 4 - E) sin /. 


* The values of b, B, c, C in the expressions for a planet’s or comet’s helio¬ 
centric coordinates referred to the equator— 

x — r sin a sin (A' 4 y), 

y — r sin b sin (B' + r), 

z — r sin c? sin (O' 4 - y), 


may be found, in a similar manner, by means of the formulse 


sin a sin E = + sin i cos & 
sin a cos E = + cos i 

cos a= 4-sin? sin £3 
sin a sin A <= 4- cos S3 
sin a cos A = — sin £3 cos i 
A'. = A + « 


sin b sin (B — A)= — cos (E+ e) 
sin b cos (B—A) = 4 - sin (E 4 - e) cos a 
cos b — — sin (E 4 e) sin a 
sin c sin (C — A) = — sin (E + e) 
sine cos (C —A) = — cos (E + e) cos a 
cos + cos (E+ e) sin a 


so that the B' 4 -y and C'+ v in the expressions for y and z are obtained by 
simply adding B — A and C — A to the values of A' 4- v. 

Are these convenient formulse for b, B, c, C not to be found anywhere in 
print? Considering that they refer to computations of such frequent occur¬ 
rence, it is not easy to assume that they have not been pointed out before. 

In case the inclination i is of moderate amount, it may be worth while, for 
greater accuracy, to determine A by computing 


sin v - tan \i . sin £3 . sin E 
which makes A = £3 4- 90 ° — v . 


If formulse for control are required, 

sin b sin B — sin S3 cos e 
sin c sin C = sin £3 sin e 
cos a 


tan (C —B) = 


cos b cos G 


may serve the purpose. 
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The positions of the satellite are then found by 

A 7 sin 5 sinp = r sin/sin (F d-^) = A.| 

A y sin s cosp = r sin g sin (Gr + u) = A.y 

Aj cos s=A + r sin h sin (H + a) = A . (I + £); 

or, if the difference between s and <d" tan's is neglected, 
rectangular coordinates by 


x" = — . • * -.sin / sin (F + u) 
A i + C 


y - 


A * 1 +£ 


sin g sin(G- + w), 


arc 1 


the 


By the substitution of a, S for A, D in the computations of the 
auxiliary angles/, F, g, G, or by computing their values, not for 
the apparent place of the planet, but for that of the satellite, the 
divisor i +£ will be got rid of, and the formulae become simpli¬ 
fied. If that course is adopted they will be 


x" = a 1 ' sin s sin p = oo" cos D sin (a — A) = a 0 vp . sin/sin (F 4- u) 

^" = sins cos p- «"sin(5 — D)— x”. tan f(a-A) sin S — a 0 vp . sin g sin(G + u) 


or, if F=F + Q, G'=G + Q, 

x" = a 0 v sin /. p gin (F' + v) 
y" = a 0 v sin g . p sin (G' + v). 


The equations of condition are 


Sx" = + v sin/ [cos (F' + v) + e cos F 7 ] . k sec p5/ 0 

— v sin /[cos (F' + v) (p cos v cos <p + e cos p + tan h <P) 4- cos F']. 

a 0 (tan p + Se) sin 5 (Q, + N) 

+ v sin / [cos (F' + v). p sin v sec 2 p - sin F']. a 0 [(tan p + Se) cos S (Q + N) 

— tan p] 

( gjv 

sin J + — J sin 5N 

( JJv 

J H—-I cos 5N-sin J] 

t" 

H- 0 

a 0 


Sy" = + v sin g [cos (G' + v) + e cos G'] . k sec p . 5/ 0 

— v sin g [cos (G ; + v) (p cos v cos (p + e cos <p + tan l 2 p; + cos G'] . 

fl rt (tnu o - 5c) sin 5 (Q + N) 

+ v sin g [cos (G' + v) . p sin v sec 2 p — sin G']. f?/(tnn p - 5c) cos 5(0, + N) 

— tan p] 

— yp [cos (Z —N) cos 5 + tanIJ sin (l — ±s + a — N) sin 5 j . 

8 /(sinJ + g).i».N 

+ vp sin {l— N). cos^. a 0 [^sin J + 5^ cos sin J] 

+ —. 5 a Q . 
a 0 

If the assumed elements are circular, and if F y = F — N, 
G v = G — 1 ST, the coordinates become 

x n = a 0 v sin/. sin (F, + 1) 
y" = a Q v sin g . sin (G, + 1) 
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and the equations of condition 
dx" = + v sin/ cos (F / 4 - 1). kM 0 

— v sin /[cos (F^^) cos 1+ cos F y ] . a Q e sin (Q + N) 

4- v sin/ [cos (F, + 1) sin l — sin F ; ] . cc 0 e cos (Q 4 - N) 

-i-ytan-^J. cos (l—N + a — N ). a 0 ^sin J 4 - sin SFT 

+ y sin (l— K) cos/. a 0 [sin^J 0 +^^cos 5FT —sin J] 

x" ~ 

+ __ . Sa Q 
a 9 

8y" = 4 * v siny cos (G y + l ). k<51 0 

— r sin g [cos (Gr, 4 - 1) cos l 4 - cos Gr,] . a Q e sin (Q, + N) 

4-J'siny [cos (Gr, 4 - ^)sin l~ sin Gr,] . a 0 -ecos(Q + N) 

— v [cos {l — FT)cosS 4 - tan^J sin (l — FT + a — N) sinS] ,« 0 sin^ 

/ S J\ 

+ y psm(l — FT)cosy . cc 0 sin[/J + ~ jcosSN —sin J] 

4- y ~~. 5& 0 
«o 

The auxiliary angles are found by the formulae 

sin/sin E= 4- sin J sin (a — N) sin y sin (G — E) = — sin (E 4 - 5) 

sin/cos E= 4 - cos J siny cos (Gr —F) = — cos(E 4 -&) cos./ 

cos / = — sin J cos (a — FT) cosy = 4 -cos (E 4 - 5) sin f 

sin/sin F = — sin (a—FT) 
sin/cos F— 4 - cos (a — N) cos J 
or F —Y by 

sin (F — N 4 - a) = tan -J. cos (a —N). sin E. 

The values of B and P used in the formulae for polar-coordinates 
are 

cos B sin P = — cos (a — FT) sin J = cos/ 
cos B cos P = 4 - cos (E 4 - 5) sin/ — cos y 
sin B = — sin (E 4 - 5) sin/ 


J 4 - — \sin5N 

w 


On the Orbit of O 2 400. By J. E. Gore. 


Some measures of this close binary star, made in 1885 with 
the 23-inch Refractor of the Princeton (U.S.A.) Observatory, and 
kindly sent me by Professor Young, show that the companion 
has described about i9o°of its apparent ellipse since its discovery 
by 0 . Struve in 1844. 

I have computed the orbit, and find the following provisional 
elements :— 


Elements of 02 400. 

P= 170*37 years £=146° 20' 

T = 1882*09 a = 43 0 30' 

^ = 0*669 tz = o / '*59 

7=36° 58' y. = — 2°*i 13 

The following table shows a comparison between these ele¬ 
ments, and the observations used in the calculation of the orbit:— 
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